Structural anomalies, spin transitions and charge disproportionation in LnCoO^ 
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The diamagnetic-paramagnetic and insulator-metal transitions in LnCoOs perovskites {Ln = La, 
Y, rare earths) are reinterpreted and modeled as a two-level excitation process. In distinction to 
previous models, the present approach can be characterized as a LS-HS-IS (low-high-intermediate 
spin) scenario. The first level is the local excitation of HS Co''+ species in the LS ground state. The 
second excitation is based on the interatomic electron transfer between the LS/HS pairs, leading 
finally to a stabilization of the metallic phase based on IS Co^^. The model parameters have been 
quantified for Ln = La, Pr and Nd samples using the powder neutron diffraction on the thermal 
expansion of Co-O bonds, that is associated with the two successive spin transitions. The same 
model is applied to interpret the magnetic susceptibility of LaCoOs and YC0O3 . 
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The perovskite cobaltites LnCoOa {Ln =- La, Y, rare- 
earth) systems show an insulating ground state based on 
the diamagnetic low spin state of trivalent cobalt (LS, 
5 = 0). With increasing temperature the systems 
undergo two magnetic transitions connected with excita- 
tions cither to the intermediate spin (IS, t\gG^g, S' = 1) or 
to the high spin state (HS, i^^e^, S = 2). 

Both transitions can be easily detected in the magnetic 
susceptibility in the case of parent compound LaCoOs. 
The first transition starts at about 50 K and saturated 
at about 150 K where a new phase is stabilized with in- 
sulating character and Curie- Weiss susceptibility. The 
saturated paramagnetic phase is manifested by a linear 
part of the inverse susceptibility over the region 150 - 
350 K that gives ^e// ~ 3.2 hb and AFM interactions 
characterized by 6* '-^ —200 K. Recent experiments at low- 
temperatures, including electron spin resonance (ESR) 
[H, inelastic neutron scattering and X-ray magnetic 
circular dichroism (XMCD) [H, identify the first transi- 
tion as an excitation to HS state. The LS/HS scenario is 
also supported by band structure calculations [1, Q that 
reveal three possible ground states for LaCoOa in the 
dependence of lattice volume and geometry: LS phase, 
LS/HS(1:1) ordered phase and homogeneous IS phase. 
The calculation evidence that the HS state is only stable 
if 6 nearest Co neighbors are in LS state, whereas sta- 
bilization of IS state is promoted by IS neighbors. HS 
states prefers antiferromagnetic interactions, whereas IS 
states prefers ferromagnetically aligned neighbors 01 . 

The second transition of an insulator-metal (I-M) kind 
starts at about 500 K and is accompanied with another 
increase of magnetic susceptibility and decrease of elec- 
tric resistivity to a metallic value p ~ 1 mficm. Above 
the I-M transition, another region with quasilinear de- 
pendence of inverse susceptibility begins at about 600 K 
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and spreads up to 1000 - 1100 K. There is a clear change 
of slope with respect to region 150 - 350 K that is tradi- 
tionally interpreted as an increase of effective moments 
of Co'^^ while the strength of AFM interactions is re- 
tained (see 6.17. Ref. 0)- We infer, however, that the 
observed susceptibility is a combined effect of a change 
of AFM interactions towards FM ones and onset of tem- 
perature independent Pauli paramagnetism, while the ef- 
fective moments remain approximately the same. 

This gives a possibility for a new interpretation of 
the spin transitions using the LS-HS-IS scenario: The 
diamagnetic-paramagnetic transition in LaCo03 origi- 
nates in thermal population of HS states in LS matrix at 
about 50 K. Each Co^+ thus appears in a mixture of these 
two states with increasing probability of the HS one that 
saturates finally in the 1:1 ratio at about 150 K because 
of strong HS/LS nearest neighbor correlations. The sec- 
ond transition (I-M) starts above 350 - 400 K and should 
be related to a charge disproportionation 0, Q , in which 
t2g electron transfer between LS/HS Co'^+ neighbors, LS 
Co3+ t%el + HS Co3+ t\^el ^ LS Go^+ tl^e^ + HS 
Co^"*" ^2#^g, with charge transfer gap ^ 0.7 eV (see Fig. 5 
in Ref. @), is an initial step. With increasing temper- 
ature, the probability of the electron exchange increases 
and new IS Co'^+ states are gradually populated due to 
recombination reaction LS Co''+ t2g^°g + HS Co^+ ^2g'^'g 
^ 2 IS Co'^+ *2g6g. This opens an Cg conduction channel 
at about 500 K and leads finally to a formation of metallic 
phase of IS character that coexists with residual regions 
of LS/HS kind with fast dynamics. The optical conduc- 
tivity data reveal indeed that this transition is associated 
with a collapse of the Cg charge transfer gap (originally 
^ 2 cV) and the effect of high temperature is thus anal- 
ogous to the effects of heavy doping in Lai_xSrxCo03 
compounds @. The transitions temperatures for other 
LnCoOs compounds are increasing with decreasing Ln 
radius [lOj. In particular, the onset of magnetic transi- 
tion in YC0O3 is shifted to 450 - 800 K and it develops 
concurrently with the I-M transition without reaching 
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FIG. 1: Experimental (symbols) and fitted (lines) thermal 
expansion of Co-0 bond for LaCoOs, PrCoOs and NdCoOa. 
The three contributions a'"", a™''^ and a™"^ of the fitted 
thermal expansion are also displayed. 



the saturated LS/HS (1:1) phase. 

In the present paper we apply the LS-HS-IS model 
to interpret the thermal expansion of Co-0 bonds of 
LnCoOs pcrovskitcs {Ln = La, Pr and Nd) and quantify 
its parameters. The transitions are treated in a rather 
formal way as a two-level excitation from LS/LS pairs 
to LS/HS and IS/IS pairs. Since paramagnetic suscepti- 
bility is not sensitive to cobalt moments in systems with 
magnetic rare earths, the present analysis is primarily 
based on observed anomalies in thermal expansion that 
accompany both transitions due to the increasing ionic 
radius of Co'^+ with increasing spin state [13, [HI, [l2, [i3| • 

The thermal expansion data have been obtained by 
neutron diffraction experiments, performed at low tem- 
perature in LLB (Saclay, France) on the G41 diffractomc- 
ter using a wavelength 2.422 A and at high tempera- 
ture in ILL (Grenoble, France) on the D20 diffractome- 
tcr using a wavelength 1.361 A. The observed neutron 
diffraction patterns were analyzed by a Rietveld method 
with the help of the FULLPROF program (Version 3.30 
- Jun2005-LLB JRC). The magnetic susceptibility was 
measured on a SQUID magnetometer in the range up to 
400 K using DC fields 10 kOc. The high temperature 
data up to 900 K were obtained using a compensated 
pendulum system MANICS in a field of 10 - 19 kOc. 

The thermal expansion a of LnCoO^ is a sum of three 
contributions: normal lattice expansion q/''^** and anoma- 
lous expansions due to the two spin transitions a™"^. 

(1) Lattice term a'"** is a weighted sum by popula- 
tions Pi{T) of thermal expansions of LnCoOa in individ- 
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FIG. 2: Upper: Population of ground and excited states. Mid- 
dle: Evolution of Weiss 6 for excited states. Lower: Molar 
susceptibility for LaCoOs and YC0O3. The solid line is fit 
based on our model (see Eq. [7|. 

ual magnetic states 
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a''^"(T) = ^ft(r) a^""(r) (1) 

i=0 

where i = corresponds to ground state LS pairs, 
I = 1 is the first excited state LS/HS pairs and i = 2 
means the second excited state IS pairs. Lattice ther- 
mal expansion increases with temperature and satu- 
rates around Debye temperature 9d- In our analysis 
we have used the saturated values ap"" = 14xl0~®K~^, 
Q,/att _ ^latt „ 18xl0~^K~^, and Debye temperatures 
Od ^ 300 K for LaCoOg and Od 700 K for the other 
LnCoOs. 

(2) The anomalous contribution a™"^ can be expressed 

as 

c^r'{T)^^jp^dr^{T) (2) 

where d™"^ is related to the difference between the size 
of ground and excited spin state and their dependence 
on temperature 

^raag^rj.^ ^ jrnag ^ a^"" (T) - (T) dT (3) 

Jo 

In our fit we have found (i™'^^(0) = 0.2% and 
d^°ff(0) = 0.7% for LaCoOa and d'^'^^iO) = 0.3% and 
4"°9(0) = 1.0% for the other LnCoOg. 

Population of excited states pi {i =1, 2) is calculated 
for each temperature point by solving the set of 2 equa- 
tions 

^ = e-^-/^ (4) 

i-E-=iP,(r) 

where Ei is the energy difference between the ground and 
excited states in units of T, and Vi is degeneracy of the 
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TABLE I: Ionic radii rLn for 9-fold coordination and the pa- 
rameters of the thermal expansion fit for Ln — La, Pr and 
Nd (see Eq. [5] and [6] for description of the parameters) and of 
the susceptibihty fit for Ln — La and Y. 
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LaCoOs 


1.216 


160 
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3.00 


2 600 


2 500 


PrCoOa 


1.179 
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370 


3.11 


2 700 


2 600 


NdCoOs 


1.163 
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3.17 


2 800 


2 650 


YC0O3 


1.075 


2 875 


860 


3.50 


3 600 


3 300 



excited state. The ground state has j/q = 1- For the first 
excited state i^i = 3 can be anticipated for spin-orbit 
spfit state of HS 00"^+ that is isolated in diamagnetic 
matrix while the spin degeneracy z^i = 25" + 1 = 5 is 
more appropriate for concentrated LS/HS pairs. For the 
second excitation 1^2 = 21^1 was used. 

The energy difference for the first excited state Ei{T) 
depends on the structure changes induced by tempera- 
ture Q or by pressure [31 ■ In our analysis the depen- 
dence of energy Ei{T) on temperature was arbitrarily 
fitted by a power function [lo| 



(5) 



where E° is the energy splitting at T = K, To is the 
temperature where Ei{To) = and n is a fitting param- 
eter that describes the curvature. 

Since the stability of IS state is promoted by neighbors 
of the same kind, the energy E2{p) was set to depend on 
the concentration p2 by the equation 0, Q 



E2{p) ^ E,{T) + E° - EP p'^ 



1/3 



(6) 



where is the energy splitting for p = and E2 de- 
scribes the dependence on p2 ■ 

The experimental and calculated thermal expansion of 
Co-0 bond for PrCoOa and NdCoOa, completed with 



data on LaCoOs taken from Ref. are displayed in 

Fig.m 

The susceptibility was fitted by an equation 



X{T) 



2 2 2 



(7) 



where Na is the Avogadro number, /i^ the Bohr mag- 
neton, Ub the Boltzmann constant and fii is the effec- 
tive moment of corresponding excited state. The Weiss 
temperature 9 was set to depend on the population of 
excited states 9i{T) = 6°pi{T). In our fit we have found 
ei = -200 K, 6if = +200 K, ^1 = /i2 = 3.25 for LaCoOg 
and 2.93 for YC0O3. 

The temperature dependence of pi and 9i for LaCoOs 
and the experimental and calculated molar susceptibility 
for LaCoOs and YC0O3 are shown in Fig.IH The param- 
eters of the thermal expansion fit (La, Pr, Nd) and of the 
susceptibility fit (La, Y) are summarized in Table HI 

In conclusion, the diamagnetic-paramagnetic and 
insulator-metal transitions in LnCoOs can be success- 
fully explained within a model of two-level excitation. 
The first level is the excitation of HS 00"^+ species in 
the LS matrix. The steep temperature decrease of the 
excitation energy close to the LS-HS crossover {To) is a 
signature that the local excitations lead finally to a global 
phase transition. The second excitation is based on the 
interatomic electron transfer and stabilization of IS Co'^^ 
states. The present interpretation can be characterized 
as a LS-HS-IS scenario that is distinct from LS-IS-HS (or 
LS-LS/HS-IS/HS) models previously used. 
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